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Magnetic excitations in ferromagnetic systems with a noncollinear ground state magnetization 
experience a fictitious magnetic field due to the equilibrium magnetic texture. Here, we investigate 
how such fictitious fields lead to thermal Hall effects in two-dimensional insulating magnets in which 
the magnetic texture is caused by spin-orbit interaction. We find that, besides the well-known 
geometric texture contribution to the fictitious magnetic field in such systems, there exists also an 
equally important but often neglected contribution due to the original spin-orbit term in the free 
energy. We consider the different possible ground states in the phase diagram of a two-dimensional 
ferromagnet with spin-orbit interaction: The spiral state and the skyrmion lattice, and find that 
thermal Hall effects can occur in certain domain walls as well as the skyrmion lattice. 
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I. INTRODUCTION 

Traditionally, the role of information carrier in spin- 
and electronic devices is taken by respectively the spin 
or the charge of the conduction electrons in the system. 
In recent years, however, there has been an increasing 
awareness that spin excitations in insulating magnets (ei- 
ther magnons or spinons) may be better suited for this 
task. The reason behind this is that these excitations 
are not subject to Joule heating. Therefore, the energy 
associated with the transport of a single unit of infor- 
mation carried by a magnon (or spinon) current can be 
much lower in such insulating systems^ The quasiparti- 
cle bosonic nature of the magnons, furthermore, allows 
in principle to essentially eliminate at low temperatures 
losses due to scattering or contact impedances^ 

Creation of a magnon current has been shown to be 
possible by means of the spin Seebeck effect^ the spin 
Hall effect*^ and with high spatial accuracy by means of 
laser-controlled local temperature gradients^ The result- 
ing spin current can be measured utilizing the inverse 
spin Hal effect It has been shown that the magnon 
current can propagate over distances of several centime- 
ters in YIG4 It has recently been shown that it is the- 
oretically possible to implement the analogs of different 
electronic components using insulating magnets^^ 

Hall effects for magnon currents are of interest both 
from a fundamental point of view as well as from the 
point of view of applications. Even though the physi- 
cal magnetic field does not directly couple to the orbital 
motion of neutral magnons, certain kinds of spin-orbit 
interactions can lead to Hall phenomena similar to that 
of a charged particle in a magnetic field. Mechanisms 
that have been shown to give rise to nonzero Hall con- 
ductances in certain insulating magnets include coupling 
of spin chirality to a magnetic field 8 and DM interaction. 9 
Of interest for applications is the fact that Hall effects in 
insulating magnets allow one to control the magnon spin 
current. 

Recently, Katsura et al. predicted 8 a nonzero ther- 



mal Hall effect for the Heisenberg model on the Kagome 
lattice. The effect originates from the fact that the cou- 
pling of spin chirality to an applied magnetic field acts as 
a fictitious magnetic flux for the magnons in the specific 
case of the Kagome lattice. Later, Onose et al. mea- 
sured" the thermal Hall effect in the pyrochlore ferro- 
magnet LU2V2O7. In this experiment, the combination 
of DM interaction and the pyrochlore structure leads to 
the finite thermal Hall conductivity. 

In those previous studies, the thermal Hall effect was 
considered using a quantum mechanical lattice model 
as starting point. The symmetry of the underlying lat- 
tice played a crucial role. We take a different approach 
and consider insulating ferromagnets with a noncollinear 
ground state magnetic texture, which we model using 
a phenomenological description. It is well known that 
the effect of the presence of a noncollinear ground state 
on the elementary excitations in a ferromagnet can be 
captured by introducing a fictitious electromagnetic po- 
tential in the equation of motion for the magnonsi 10 i n 
As we will show, there are two contributions to this elec- 
tromagnetic potential. The first one only depends on the 
equilibrium magnetic texture, the second depends on the 
form of the free energy that gives rise to the noncollinear 
ground state in the first place, i.e. the contribution to 
the free energy due to spin-orbit interaction. The latter is 
often unjustifiably neglected, but we will show here that 
it has important consequences. Indeed, since both con- 
tributions are determined by the spin-orbit interaction, 
they will be of similar magnitude. 

It has been shown that the fictitious electromag- 
netic potential couples the motion of magnetic texture 
and that of heat currents. = This coupling reveals itself 
through local cooling by magnetic texture dynamics 12 
and thermally induced motion of magnetic textures such 
as domain wallsj 13 ' 14 

This work is organized as follows: In Sec. [II] we intro- 
duce our system and derive the fictitious electromagnetic 
vector potential that acts on the magnons, which turns 
out to include an often-overlooked contribution. In Sec. 
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FIG. 1. Pictorial representation of the thermal spin Hall ef- 
fect. A temperature difference AT applied to a sample leads 
to a finite heat current. Since the heat current is carried 
by the magnons in the system, the fictitious magnetic field 
that magnons experience due to a non-trivial magnetic ground 
state will lead to a finite thermal Hall conductivity. 



IIIII we derive the relevant ground state properties of the 
different ground states in the phase diagram of the insu- 
lating ferromagnet with nonzero Dzyaloshinskii-Moriya 
interaction. In Sec. IIVI we calculate the band structure 
of one of the ground states, the triangular skyrmion lat- 
tice, and calculate its thermal Hall conductivity. 



II. MAGNONS IN THE PRESENCE OF 
MAGNETIC TEXTURE 

We consider a two-dimensional non-itinerant ferromag- 
net in the x-y plane with spatially varying and time de- 
pendent spin density sm(r, t). The spin density is related 
to the magnetization M(r,t) as sm(r,f) = M(r, i)/7, 
where 7 is the gyromagnetic ratio (7 < for electrons). 
The magnitude s of the spin density is assumed to be 
constant, and m(r, t) is a unit vector. The system is 
described by the Lagrangia n 10 ' 11 



d 2 r [D(m) • rh — J-(m, djin)] 



(1) 



Here D = sh (n x m) /(l + m • n) is the vector potential 
corresponding to Wess-Zumino action with an arbitrary 
n pointing along the Dirac string, J-(m, djm) is the mag- 
netic free energy density of the system, which we assume 
to be of the form (double indices are summed over) 

F{m,djia) — — (djm.) 2 — M s m-H+sFr(m, djtn). (2) 

Here J is the strength of the exchange interaction, M s — 
7s is the saturation magnetization, H the external mag- 
netic field (which we will always assume to be in the z di- 
rection), and J-r(m, c^rri) describes terms due to broken 
symmetries. For isotropic ferromagnets in the exchange 



approximation, the leading order terms in the free en- 
ergy are quadratic in the texture [first term in Eq. (JTJ)]. 
Breaking inversion symmetry by spin-orbit interactions, 
while still retaining isotropy in the x-y plane, allows to 
construct terms that are first order in texture. These 
terms are given by 

Jr(m, djin) = r R m z V • m + Tdmhi • (V x m) . (3) 

Here, we defined V = d x x + d y y. The first term is due 
to structural inversion symmetry breaking and hence is 
anisotropic in the z direction. Such terms occur in sys- 
tems with finite Rashba spin-orbit interaction^ or on the 
surface of a topological insulator^ The second term de- 
scribes Dzyaloshinskii-Moriya interaction^ which origi- 
nates from the breaking of bulk inversion symmetry and 
is therefore isotropic. We note that the two terms in Eq. 
Q are equivalent (up to an irrelevant boundary term) 
under a simple rotation around the z axis in spin space. 
Since all other terms in Eq. ([!} are invariant under such 
rotations, we can always absorb the term proportional 
to Tr in the term proportional to Tdm- We will there- 
fore put Tr to zero in the remainder of this work. For 
simplicity, we have ignored a term —nm 2 z that would de- 
scribe easy axis anisotropy, and a term —M s m • H m /2, 
where H m describes the magnetic stray field, in Eq. ©. 
Substitution of Eq. ([T]) in the Euler-Lagrange equation 
leads to the Landau-Lifshitz equation 



shm — m x 5 m F(m, djin) = 0, 



(4) 



where F(m, djm) is the total free energy of the system. 
We split the magnetization m in a static equilibrium 
magnetization mo and small fast oscillations 5m (spin 
waves) around the equilibrium magnetization. To lowest 
order in 5m the two are orthogonal. In a textured magnet 
mo = mo(r), which makes finding the elementary exci- 
tations a nontrivial task. To circumvent this issue we in- 
troduce a coordinate transformation m'(r) = i?(r)m(r), 
where R(r) is such that the new equilibrium magnetiza- 
tion m is constant and parallel to the z axis. In this 
coordinate frame the spin waves are in the x-y plane. 

The 3x3 matrix R describes a rotation over an an- 
gle 7T around the axis defined by the unit vector n = 
[z + mo] / [2 cos {0/2)]. Here, 6 is the polar angle of m . 
Using Rodriques' rotation formula, we findfl= 2nn T -i. 

The effect of the transformation to the new coordi- 
nate system is that we have to use the covariant form of 
the differential operators, 9 M — > (9 M + A^), with A M = 
R (d^R), in the Landau-Lifshitz equation. The sub- 
script /1 describes both time (fi — 0) and space (fx = 1,2) 
coordinates. 

In the new coordinate system, the Landau-Lifshitz 
equation for the free action Eq. (J5]) becomes 



ihdtm 



J (V/i + A) +ip m+ 



(5) 



Here, m± = (5m! x ± i5m' y )/\/2 describe circular spin 
waves in the rotated frame. Furthermore, ip — mo-H/s+ 
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h\Er 1 (d t R)\ 12, and the components of the vector poten- 
tial A are given by Aj = Aj. 12- The skew-symmetric 
matrix A is here defined as Aj = R{dj — QIj)R. In the 
latter equation we defined C = ^t>m/Ji and the matrices 



(6) 



We see that the components Aj of the fictitious mag- 
netic vector potential consist of two contributions. The 
first comes from the exchange interaction in the presence 
of magnetic texture; the second (texture-independent) 
part originates from the DM interaction term in the free 
energy. The latter contribution is often unjustifiably ne- 
glected, but we will show here that it has important con- 
sequences. Indeed, typically both contributions will be of 
the same order of magnitude. This is because the mag- 
netic texture itself is caused by the DM interaction, and 
will therefore manifest itself on lengthscales J/r^M- 

We can requantize Eq. ([5]) by introducing the bosonic 
creation operator b* oc m_. This requantization works 
since m' + and m'_ satisfy approximate bosonic commu- 
tation relations in the limit of small deviations from equi- 
librium. After requantization, Eq. ([5]) can be interpreted 
as the von Neumann equation belonging to the Hamilto- 
nian 



H = / d 2 r^ 



J(V/i + A) 2 +if 



(7) 



Therefore, the elementary excitations of the system be- 
have as noninteracting bosonic quasiparticles which in- 
teract with the collective smoothly-varying background 
that is treated classically by the inclusion of a fictitious 
magnetic vector potential A and electric potential if. 

In the derivation of Eq. ([5]) we have assumed that 
the length of a typical wave packet is much smaller than 
the spatial extension over which the magnetic texture 
varies. We will refer to this as the adiabatic approxima- 
tion^ Using this assumption, we have neglected terms 
that are higher order in texture in Eq. ([5]). Such terms, 
which become important at lower wave vectors, con- 
tribute to the fictitious electric potential ip in Eq ([5]) 
a term - J[i 2 , 13 + A\, 2Z + i 2 43 + i 2 , 23 ]/2, which is 
quadratic in magnetic texture. Furthermore, at lower 
wave vectors one needs to add to the right-hand side of 
Eq. © aterm J[(A x , 13 +iA x a3) 2 + (Ay,i3+iAy,23) 2 }rn^, 
which describes the finite ellipticity of the magnons. 



III. TEXTURED GROUND STATES 

In this section we will present the different possible 
ground states for systems with free energy given by Eq. 
© (with Tr = 0), as a function of external magnetic field 
H = Hz. We will also present the fictitious magnetic 
vector potentials that find their origin in these textured 



ground states. It has been showniS that as the mag- 
netic field H increases from zero, the ground state of the 
two-dimensional fcrromagnet with spin-orbit interaction 
changes from a spiral state for H < H c \, to a skyrmion 
lattice state for magnetic fields H c \ < H < H C 2, and 
finally to the ferromagnetic ground state for H > H C 2- 
Both critical fields H c \ and H c2 are typically of the order 



Fq M /J (see Ref. [20j). This last observation, in combina- 
tion with the adiabatic assumption and the fact that the 
equilibrium magnetization is time- independent, allows us 
to neglect the fictitious electric potential ip in Eq. ([7|). 
Since the ferromagnetic ground state has no magnetic 
texture, it is of no interest for our purposes. In this sec- 
tion we will therefore derive the properties of the spiral- 
and skyrmion lattice ground state. 



A. Spiral state 

Following Ref. [l9| we will derive the properties of the 
spiral ground state in a 2-dimensional ferromagnet with 
DM interaction. We write mo(r) in the following form 

mo(r) = cos£sin6>x + sin^sin^y + cos6*z. (8) 

For the spiral state, 9 = 9(y) and £ is a constant. With 
these constraints, the free energy becomes a functional 
that depends only on 9(y) and d y 9(y). Minimizing this 
functional with respect to 9(y) gives the following differ- 
ential equation 



d 2 y 9 + asm9 = 0, 



(9) 



where we defined a = —jH/ J. Eq. is the equation 
of motion for the mathematical pendulum. The general 
solution is given in implicit form by the expression 



Hv) d9 



1 



1 - to 2 sin 2 9/2 



(10) 



where to = 4a/ (2a + C) and f3 = \j2a + C. Alterna- 
tively, we can write 9(y) — 2<j)(f3y/2, to), where </)(it,m) 
is the amplitude of the Jacobi elliptic function. The con- 
stant C is the first constant of integration. To determine 
it, we use the fact that 9(y) is a periodic function with 
period y^. By integrating the inverse of the first integral 
d y 9 of Eq. © over one period we can determine yo as 



yo 



2tv 



d9- 



1 



y/2a cos 9 + C 



(11) 



To fix C, we minimize the average free energy 
(l/j/o) Jq° d y 9), which leads to the following implicit 
expression for C 



p2tt 

/ d9y/2acos9 + C = 2ir(. 
Jo 



(12) 



The minimization of the average free energy also fixes 
cos£ = 1. From this we see that the ground state is 
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FIG. 2. Fictitious magnetic field due to the spiral ground 
state. Parameters are C = 70 jtm"\ J/(fese 2 ) = 63 K, and 
the interatomic spacing is taken to be e = 4.5 A. (See Ref. l2(il ) 
To make the connection to electromagnetism, we note that a 
fictitious field B z = 27r£/yo(0) ~ 5 ■ 10 15 m -2 acting on a spin 



wave gives rise to the same magnetic length as a 
magnetic field acting on a free electron. 
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a spiral state whose structure locally resembles a Bloch 
domain wall, as is expected for the DM interaction. — We 
also note that in the case of zero magnetic field (a = 0) 
the spiral state is described by a simple sinusoid with 
period yo = 2tt/(, whereas for finite magnetic field the 
mirror symmetry with respect to the z plane is broken. 
Eq. (|12[) also puts a constraint on the maximum value of 
H for which the spiral state is stable. 

Some general observations can be made with regards to 
the fictitious magnetic vector potential due to the spiral 
ground state. For the general ground state Eq. ([8|) with 
9 = 9(y) and £ constant, the fictitious vector potential is 
A(y) — ( sin 9(y)ic. This potential is caused solely by the 
DM contribution to A, the geometric texture contribu- 
tion is zero everywhere. The z component of the fictitious 
magnetic field that the magnons experience is given by 
B z {y) = V x A| z = — C,d y 9{y) cos9(y). It is easily seen 
that the total fictitious magnetic flux through one pe- 
riod of the spiral (B z ) = dyB z (y) = 0. The fictitious 
magnetic field B z (y) has been plotted in Fig. [2]for differ- 
ent values of the applied magnetic field H. Transport in 
the presence of a magnetic field that is nonhomogeneous 
in one direction and has zero average has been studied 
extensively (see Ref. [2l| for a recent review). It is well 
known that these systems do not display a finite Hall 
conductivity. However, such magnetic fields have been 
predicted to influence the longitudinal conductance, due 
to the presence of localized snake orbits at energies that 
are low compared to the cyclotron frequency associated 
with the amplitude of the magnetic field . 22 ' 23 From our 
analysis it is also seen that one dimensional textures can 
give rise to a nonzero average fictitious magnetic flux for 
certain domain walls, since these consist of only half a pe- 
riod of the spiral. Hence, such domain walls will display 
the thermal Hall effect. 



B. Skyrmion lattice 

For magnetic fields H c \ < H < H c2 the ground state 
of the two-dimenional ferromagnet with DM interaction 
is a skyrmion latticed This triangular lattice has ba- 
sis vectors ai = ax. and a.2 = (a/2)x + (a-\/3/2)y, and 
contains skyrmions with radius R. The size of a single 
unit cell is (v / 3/2)a 2 , where a = 2R. The magnetization 
m (r) of a single skyrmion of radius R is parametrized 
in polar coordinates (p, <fi) by Eq. © with 9 = 9(p) and 
boundary conditions 



(9(0) = 7T and 6(R) = 0. 



(13) 



Furthermore, £ = N<f> — ir/2, where N is the charge of 
the skyrmion. We will assume N = 1 throughout. The 
magnetization profile can in principle be determined nu- 
merically by minimizing the free energy with the afore- 
mentioned boundary conditions. However, for simplicity 
we will assume a linear dependence 9(p) = 7r(l — p/R) 
for our principal analysis of the texture-induced thermal 
Hall effect. 

In polar coordinates the fictitious magnetic vector po- 
tential A(r) due to a single skyrmion is given by (here 
< p < R) 



A(r) 



cos 9(p) — 1 



(cos9(p) 



(14) 



The z component of the fictitious magnetic field for this 
vector potential is given by B z (p) = p^ 1 d p {pA t j } ). It fol- 
lows that the total flux through a unit cell is (B z ) = 
2tt J^ 1 dppB z = 4-7T. This means that each unit cell con- 
tains two magnetic flux quanta. The nonzero average 
flux is caused by the texture contribution to A(r), the 
DM-intcraction contribution averages to zero. From the 
fact that the average magnetic flux is nonzero, it follows 
that the skyrmion lattice has a nonzero Hall conductiv- 
ity. One might then be inclined to take the average values 
of the fictitious magnetic field and ignore the spatial de- 
pendence when calculating the thermal Hall conductivity 
of the skyrmion lattice. However, we will show shortly 
that the spatial variation of the fictitious magnetic field 
is substantial, so that we should take both contributions 
into account in our analysis. 

Let us consider for simplicity the situation in which 
R = 7r/£. In that case B z (p) = cos 9(p). The spa- 
tial variation is therefore large enough that the fictitious 
field switches from a negative minimum at p = to a 
positive maximum at p = R. Such large variations have 
been shown to have a significant influence on the band 
structure of magnetic lattices^ 

For what follows, it will be convenient to formally 
split the fictitious magnetic vector potential in two parts, 
A(r) = Ao(r) + A'(r), where Ao(r) describes the contri- 
bution from the homogeneous nonzero average fictitious 
magnetic flux, and A'(r) the periodic contribution with 
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zero average (we work in the Landau gauge) 
A (r) = -B yx 

A'(r) = J2 [Mr, r?)* + A y (r, r?)y] e<^+^)- r .(15) 

Here, i?o = 87r/(\/3a 2 ) is the average fictitious mag- 
netic field, and = (2n/a)(it — y/y/3) and k 2 = 
(27r/a)(2/v / 3)y are the basis vectors of the reciprocal lat- 
tice, such that the periodic part of the fictitious vector 
potential satisfies A'(r + ai) = A'(r + a2) = A'(r). Such 
spatially varying magnetic fields are known to give rise 
to finite Hall effects, even in the absence of a nonzero 
average.— 



IV. THERMAL HALL CONDUCTIVITY OF 
THE SKYRMION LATTICE 

Since the skyrmion lattice can be described by a free 
bosonic Hamiltonian with a spatially varying fictitious 
magnetic field with on average two magnetic flux quanta 
per unit cell and the same symmetry as the skyrmion 
lattice, the eigenstates of the skyrmion lattice are mag- 
netic Bloch states. In Sec. IIV Al we will determine the 
excitation spectrum and explicit form of these states. In 
Sec. IIV Bl we will show how the thermal Hall conduc- 
tivity of the skyrmion lattice is determined by the Berry 
curvature of these magnetic Bloch states. 



A. Diagonalization 

To find the elementary excitations of the skyrmion lat- 
tice, we need to diagonalize the Hamiltonian H in Eq. 
([7]) with the fictitious magnetic vector potential given in 
Eq. (fl4j) . We do this by numerically diagonalizing the 
matrix that results from rewriting H in the basis of the 
Landau levels that describe excitations with the appro- 
priate symmetry in the presence of the fictitious magnetic 
vector pot ential Ao(r) only. Our derivation follows that 
of Ref. [2J, with the difference that we consider the case 
with two instead of one flux quantum per unit cell. 

Eigenstates for the free system with only a homoge- 
neous magnetic field BqZ and without any underlying 
symmetries are given by 

AT 1 _i 

V>„ fcx (r) = -^e-* k * x v n (B%y + B ~ 2 k x ), (16) 

where N n = (Si)* and <p n {x) = ' 2 H n {x), 

with H n (x) the n-th Hermite polynomial. The corre- 
sponding energies are E n = 2JB (n + 1/2). To account 
for the presence of the triangular lattice, and the fact 
that every unit cell contains two flux quanta, we need to 
find the most general linear combination of eigenstates 
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(0,0) (1,1/2) (1,0) (0,0) 

FIG. 3. Band structure of the skyrmion lattice with parame- 
ters R = 45 nm, ( = 70 fim' 1 , and 2JB /k B « 50 mK. The 
labels on the horizontal axis denote (fci,fe2), with the wave 
vectors normalized to 2n ja. 

that satisfies 

M ai Vwk(r) = e ikia ip nmk (r), 
Ma 2 ?/Wk(r) - e jfc2 > nmk (r). (17) 

Here, k\ and fc 2 are defined such that (27r/a)k = fciki + 
£^2- Furthermore, k is restricted to lie within the first 
Brillouin zone. We will discuss the origin of the quan- 
tum number m later. We have to work with magnetic 
translation operators M ai 2 since the canonical momen- 
tum is no longer a good quantum number in the pres- 
ence of the vector potential Ao(r). These magnetic 
translation operators are defined as M ai = T ai and 
M a2 = cxp[— i(47r/a)a;]T a2 , where T ai 2 are the usual 
translation operators. The appropriate eigenstates are 
then given by 

oo 

Vw(r) - ]T (-l)(<+¥)0+f-i) e -*(*+f X^-feOa 

/— — OO 

xV> n ,_ fcl _ (i+¥) i2L. (18) 

The quantum number m, which in our case can take val- 
ues or 1, accounts for the fact that in the presence of a 
natural number p of flux quanta per unit cell each mag- 
netic band will split up in p subbands. These subbands 
are degenerate for a constant magnetic field, but will in 
general split for a spatially varying magnetic field, as we 
will see later. The set of wave functions defined in Eq. 
(1181) constitutes a complete orthonormal basis with trian- 
gular symmetry. The eigenfunctions are chosen in such a 
way that perturbations in the fictitious magnetic vector 
potential that are periodic in the triangular lattice are 
diagonal in the momenta k\ and fe. 

We are now in a position to calculate the matrix el- 
ements of H with respect to the basis defined by the 
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eigenstates in Eq. (fl"8|) . We rewrite H = H + Hi + H 2 . 
where the subscript denotes the order in which A'(r) oc- 
curs in the respective term. Hq is then trivially given 
by 

(n',m'\H \n,m) = 2JB (n + l/2)5 n , n ,5 m , m ,. (19) 
The matrix elements of H\ are given by 

(n',m'|-ffi|n,m)„'>„ = J6 m >- m<T y~] B(t, rf) 
i / -L jjj , 2tt/ / 

L n (z T7] ) — I - L n {Zrr}) ~ ~ L n _i (z Tr} ) 

\ Z-TTj ^TT] 

x{-l) m ^G n . n {r,r,), (20) 
and the matrix elements of H 2 by 

(n\ m'\Hi\n, m) n >> n = J8 m >- m ^ T > +T 
x [A x (t' , r)')A x (r, rj) + A y (r' , r)')Ay(r, 77)] 

T 1 .7]' ,T,r/ 

X (-l) m «+ri) Gn , n (T' +T, V '+ 7]). (21) 



We defined the function 

1/2 



.277 - 



-z Tr ,/2 Trirri/2 ir]kia/2 ig(fc 2 a+7r)/2 



Furthermore, we defined z TJ) = (27r/\/3)(T 2 — rrj + r\" 
The function L%(x) is the associated Laguerre polyno- 
mial. The first 10 subbands of the band structure of the 
skyrmion lattice with parameters 2J_Bo/&s ~ 50 mK, 
R — 45 nm, and C = 70 /1m -1 are given in Fig. [3] In 
our numerical calculation we used the fact that the cou- 
pling between two band decays super-exponentially [to 
be precise, it decays as y/ (n!/n'!)], so that only a limited 
number of bands have to be taken into account. It is 
seen that the inclusion of the spatially varying fictitious 
magnetic field has a pronounced effect, leading both to 
different splittings of the different subbands, as well as 
substantial broadening of the subbands. From Fig. [3] it 
is seen that the typical level splitting between magnetic 
subbands is 50 mK, which sets the temperature scale on 
which the system is in the quantum Hall regime. Sys- 
tems with larger ratio r^ M / J will display quantum Hall 
behavior at higher temperatures. We note that finite 
Gilbert damping a will broaden the different magnetic 
subbands by an amount (Aw/w) = 2a. Eventually this 
will destroy the visibility of individual subbands. How- 
ever, since the Gilbert damping is around a ~ 10~ 3 in a 
range of different materials, this only becomes problem- 
atic at high magnetic subbands. 

We note that within our model we do not find the ex- 
pected Goldstone modes associated with the skyrmion 
latticed We argue that this is due to our adiabatic as- 
sumption, which breaks down for the smallest wave vec- 
tors. Assuming a quadratic dispersion for the magnons, 
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FIG. 4. Berry curvature of the two highest magnetic sub- 
bands in Fig. [3] in a single Brillouin zone. The subband cor- 
responding to the top figure does not carry a net curvature, 
the bottom figure carries 2n. 



we can estimate the magnitude |k m | of the characteris- 
tic wave vector of the magnons that make up the lowest 
magnetic subband as J|k m | 2 = J Bo, which leads to a 
typical magnon wave length A m ~ a. The wave vector 
|k m | increases for higher subbands. Since the accuracy 
of our model increases with increasing wave vector, our 
description improves for higher magnetic subbands. 

In the next section we will investigate the effect of the 
finite bandwidth of the magnetic subbands on the ther- 
mal Hall conductivity of the skyrmion lattice. 



B. Thermal Hall conductivity 

It is well known 27 that the semi-classical dynamics of 
a wave packet in the basis of the magnetic Bloch states 



Unk(r) 



-ik-r 



^„k(r) is given by 
r = <9k£Vi(k) — k x fi„(k) and hk = 0. 



(23) 



We have assumed here that there are no electric fields and 
that the states u n k(r) are the eigenstates of the Hamil- 
tonian H including the fictitious magnetic vector poten- 
tial A(r). fi n (k) is the Berry curvature of the magnetic 
Bloch band. Since we consider a two-dimensional system, 
only its z component is relevant. It is given by 



n„(k) = 2Im 



"nk(r) 

dk x 



u„k(r) 



dk v 



(24) 



For the skyrmion lattice, the magnetic Bloch states are 
given by 



Unk(r) = e lkr ^ c", m , k i/wk(r). 



(25) 



n' ,m' 
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The weights c", m , k follow from the diagonalization per- 
formed in Sec. IIV Al It should be noted that every com- 
pletely filled subband carries a total Berry curvature that 
is a multiple of 27r, in accordance with the quantization 
of the TKNN invariant^ 

Matsumotu et al. have showr*2£ that the thermal Hall 
conductivity for a system described by Eq. (|23[) is given 
by 

Kxy = -JTy- 2_, c 2(P«k)^n(k). (26) 
n,k 

Here c 2 {p) = (1 + p) (log^)* - (logp) 2 - 2Li a (-p) 

describes the effect of the thermal distribution of the 
magnons, and p(e) = (exp^ 6- ^ — is the Bose dis- 
tribution function. 

Fig. 2] shows the Berry curvature for the two high- 
est magnetic subbands in Fig. [31 Together with the 
band structure given in Fig. [3] the Berry curvature com- 
pletely determines the thermal Hall conductivity, as can 
be seen from Eq. (|2l))) . The position and width of the 
magnetic subbands determine at which energies states 
become available for thermal transport, the Berry cur- 
vature determines the extend to which these states con- 
tribute to the thermal Hall conductivity. 

We have shown then that the spatially varying ficti- 
tious magnetic field gives rise to a nontrivial structure 
of the Berry curvature as well as a broadening of the 
magnetic subbands. As follows from Eq. (|21)]) . the com- 
bination of these two effects may be studied experimen- 
tally by measuring n xy as a function of temperature. As 
the temperature increases, the thermal distribution of 
the bosonic magnons broadens, which enables the higher 
bands to contribute to thermal transport. 

Alternatively, one can probe the chiral subgap edge 
states (above the bands whose total TKNN number is 
odd) directly by microwave means. The first magnetic 
subband corresponds to several GHz, which falls within 
the scope of standard microwave techniques. Excitation 
of some of the low-energy modes of the skyrmion lattice 
by microwave radiation has been performed 3 -^ and anal- 
ysed^ The magnonic reflectionless wave guide modes are 
analogous to those found in photonic crystals 3 -^ and can 
provide intriguing spintronics applications. Magnonic 
edge states have recently been proposed to exist in YIG 
without any equilibrium spin texture but with an ar- 



ray with an array of Fe-filled pillars, in which ellipticity 
of magnons due to dipolar interactions can give rise to 
magnonic bands with a nonzero TKNN invariant. 33 



V. CONCLUSIONS 

We have studied how fictitious magnetic fields, which 
are caused by textured equilibrium magnetization, lead 
to thermal Hall effects in two-dimensional insulating 
magnets in which the nontrivial equilibrium magnetiza- 
tion is caused by spin-orbit interaction. We have given a 
general expression for the fictitious magnetic vector po- 
tential, and found that it consists of two contributions: 
A geometric texture contribution and an often neglected 
contribution due to the original spin-orbit term in the 
free energy. We have shown that both contributions are 
generally of the same order of magnitude. 

We have derived the relevant properties of the two 
ground states (the spiral state and the skyrmion lat- 
tice state) in the phase diagram of a two-dimensional 
non-itinerant ferromagnet with nonzero Dzyaloshinskii- 
Moriya interaction. We have found that a system which 
has the spiral state as magnetic ground state does not 
have a finite thermal Hall conducitivity. However, we 
predicted that certain domain wall structures do display 
thermal Hall effects. 

Wc have numerically diagonalized the Hamiltonian de- 
scribing the triangular skyrmion lattice. We found that 
due to the spatially varying fictitious magnetic vector 
potential, the excitation spectrum consists of broadened 
magnetic subbands. We have calculated the Berry cur- 
vature of the magnetic subbands, and showed that the 
Berry curvature in combination with the excitation spec- 
trum completely determines the thermal Hall conductiv- 
ity of the skyrmion lattice. 
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